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Abstract 

The properties of the two-chain Hubbard Model doped away from half-filling 

are investigated. The charge gap is found to vanish, but a finite spin gap 

exists over a range of interchain hopping strength t±. In this range, there are 

modified d x 2_ y 2-like pairing correlations whose strength is correlated with the 

size of the spin gap. It is found that the pair field correlations are enhanced 

by the onsite Coulomb interaction U. 
PACS Numbers: 74.20Hi, 75.10Lp, 71.10+x 
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Studies of two chain t—J and Hubbard models provide insight into the development 
of pairing correlations in strongly coupled electron systems which can exhibit a spin-gap 
phase. For example, Lanczos calculations for a t—J model by Dagotto et al. |T] suggested 
that antiferromagnetic S — 1/2 coupled chains such as (VO)2P20y should have a spin 
gap and that if this system could be doped, it would have enhanced superconducting or 
charge-density-wave correlations in the ground state. Rice and coworkers |3],^[] discussed 
Sr 2 Cu406, which consists of layers containing weakly coupled CuO ladders ||, and based 
upon their analysis of a t—J model suggested that if this system were lightly doped, the spin 
gap phase would persist and a ground state with dominant superconducting correlations 
could be realized. 

Various approaches to the two-chain Hubbard model have also been reported |6|-[TT 



Renormalization group calculations |8|-|11]1 find evidence for a variety of phases. However, a 
number of these phases correspond to renormalization flows to strong coupling, making it 
important to have the guidance provided by independent numerical calculations. Previous 
density matrix numerical renormalization group calculations |T2[ suggested that the two 



chain Hubbard model has a spin gap at half-filling, provided the onsite Coulomb interaction 
U and the interchain hopping t± are finite. For U = 8t and t± = t, with t the one electron 
hopping along a chain, a spin gap was also found for band fillings of (n) = 0.875 and 
(n) = 0.75 electrons per site. In addition, for the doped case, the equal time pair field 
correlations were observed to decay as a power law, but with an exponent similar to that of 
the noninteracting (U = 0) system. Here we extend these two-chain Hubbard calculations 
to explore the dependence of the doped state on the strength of the interchain hopping t±. 
We find that the charge gap vanishes in the doped state, but that a spin gap exists over a 
range of t±. Furthermore, the decay of the pair field correlations depends on t±, with the 
strongest correlations appearing when the spin gap is a maximum. 
The two chain Hubbard model has a Hamiltonian 

H = -tY,( C l\a C i+lM + 4+lM C iM) ~ t± Y,( C iM C ^ + C UaCi,la) + U £ Tl iM n iM . (1) 
i.Xcr i,a i.X 



Here c iXa creates an electron of spin a on rung j and chain A, with A = 1 or 2. The 
one-electron hopping along a chain is t, the hopping between chains is t±, and the onsite 
Coulomb interaction is U. We will measure energies in units of t. 

The calculations were carried out using the density matrix renormalization group 
(DMRG) [[14]]. This technique allows the calculation of ground state energies and corre- 
lation functions in a controlled way. The details of the techniques used for two Hubbard 
chains are described in Ref. fll5]| . We studied lattices with open boundary conditions of 
various sizes up to 2 x 40 sites, keeping up to 500 states per block. The discarded weight 
of the density matrix (the truncation error) ranged from 3.7 x 10~ 5 to less than 10 -8 . We 
estimate that the maximum errors on the quantities shown in this paper are at most a few 
percent, and typically are of the order of the plotting symbol size or less. 

The charge and spin gaps are defined in terms of the ground state energy E (N-^, N^) for 
iVf spin-up and iVj_ spin-down electrons. The charge gap is given by 

A c = [E (N -l,N-l) + E (N + 1, N + 1) - 2E (N, N)]/2, (2) 

while the spin gap is 

A S = E (N + 1,N-1)- E {N, N). (3) 

These quantities are shown in Fig. [I], calculated for Hubbard ladders of different lengths L. 
In Fig. |l](a) we show Ac for U = 8, t± = 1, and several different band fillings. At half- 
filling, Ac scales to a large finite value A C (L = oo) « 4.5, whereas in the doped case, with 
(n) slightly less than 1, we find that the charge gap vanishes in the thermodynamic limit. 
For the spin gap, shown Fig. |l](b), the finite size scaling must be done carefully in order to 
determine the behavior because the finite size effects are relatively large and have a different 
character in different regimes. We fit the data using a least squares fit to polynomials in 
L _1 . At half-filling, we use a third order polynomial to fit six system sizes. At (n) = 0.875, 
fewer sizes are available with the correct filling, so we fit to a second order polynomial. The 
polynomial fits are shown as the lines in |l](b). At half-filling, Ag(oo) ~ 0.12 and is clearly 
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nonvanishing. The system is an insulator with short-range antiferromagnetic correlations. 
For the doped system, there are three distinct regimes for which we show representative 
curves with different finite size behavior. For small t±, illustrated by the t± = 0.5 curve, 
the finite size effects are large with a large linear coefficient. In the thermodynamic limit, 
the spin gap becomes small. For a smaller t± of 0.3, we find As(32) «3x 10~ 3 . Therefore, 
we believe that As(oo) vanishes for t± < 0.5. For intermediate t±, shown for t± = 1.0, 
the coefficient of the quadratic term is large and positive and the spin gap scales to a finite 
value. For large t±, shown for t± = 2.0, the finite size effects are smaller, the scaling is linear 
in 1/L and A,s(oo) vanishes. 

We plot the spin gap extrapolated to the thermodynamic limit as a function of t± in 
Fig. ||. It appears that in the doped system, a finite value of t± is required to produce an 
interchain rung exchange J± ~ At\/U which is sufficient to create a spin gap. This differs 
from the half-filled case where our earlier work suggests that a spin gap opens for all t± > 
IT^JIBfl . For t± > 1.7, as shown in Fig. |2|, the spin gap is suppressed. In a weak coupling 



picture, a two chain model leads to bonding and antibonding bands split by the coupling t±. 
For U = 0, the antibonding band is completely unoccupied at (n) = 0.875 when t± > 1.85, 
suggesting that the transition at t± ~ 1.7 is related to the band splitting. In the spin gap 
phase, the spin-spin correlations are short-range, incommensurate, and antiferromagnetic 
in nature. They give rise to a Lorentzian peak in the magnetic structure factor ||12|| . 

In order to further characterize the spin gap region, we have calculated the equal-time 
rung-rung pair field correlation function (AjAj). Here 

A J = H.r'k " 4il c l2 T ) ( 4 ) 

creates a spin singlet pair across the jth rung and A« destroys a singlet pair on the ith 
rung. For comparison with the spin gap, Fig. |2| also shows the pair field correlation function 
(AjAj) averaged over separations 8 < \i — j\ < 12. It is clear that the size of the spin gap 
and the pair field are correlated. 

Figure |3|(a) shows a log-log plot of (AjAj) versus \i — j\ with (n) = 0.875, U — 8 and 



t± = 1.5 for various length chains. As discussed earlier, the calculational procedure which 
we have implemented has open end boundary conditions. The correlation functions shown 
in Fig. ^(a) are obtained by averaging over a number of i and j for a given \i — j\. We have 
found that this procedure removes most of the finite size effects due to the open boundary 
conditions. One can see that the boundary effects are small in Fig. 0(a) by the overlap of 
the numerical results for chains of different lengths. Only when \i — j\ is large enough that 
both i and j are within a few lattice spacings of the boundaries do significant differences 
appear. The dashed curve corresponds to the \i — j\~ 2 decay of the non-interacting pair field 
correlations and one sees that for t± = 1.5, the pair field correlations of the interacting system 
are enhanced over the U = system. In Fig. 0(b), we show (AjAj) versus \i — j\ calculated 
on a 2 x 32 system for (n) = 0.875, U = 8, and various values of For small values of 
t±, the pair field correlations decay rapidly. As t± increases, the pair field correlations are 



enhanced until t± ~ 1.7, after which they are suppressed by further increasing t± |T3 . 

It is also possible to determine the internal structure of a pair. We have calculated the 
pair wave function obtained from the off-diagonal ground state expectation value 

V(iJ) = N-2(i>o\{ c rt c H ~ CjlCi^o) N - ( 5 ) 

Here {c^c^ — cj^c^) removes a singlet pair with one electron on the ith site and one on 
the jth site. Near the center of the ladder, + x,i) = 0.1 and +y,i) = —0.08 for 
U = 8, t± = 1.5, and (n) = 0.875. In the cuprate chain systems, t± and t have the same 
sign so that the relative negative sign of + x,i) and ^>(i + y, i) is physically relevant 
and corresponds to a modified d x 2_ y 2-\ike state as previously discussed P,[12|. The small 
amplitude for removing a pair from a rung, 0.1, reflects the spatial extent of the relative, 
internal coordinate of the pair wave function and the fact that the cj Acr operators are bare 
quasi-particle operators. When one takes out a factor of (0.1) 2 associated with this overlap 
of Aj and A] with the internal pair wave function, it is clear that there are significant 
pair-field-pair-field center of mass correlations present in the doped system for t± = 1.5. 
Recently, we described an intuitive picture for the spin-liquid state of two coupled chains 



16] based on a resonating valence bond (RVB) variational ansatz |17]]. This picture is also 



useful for describing the apparent destruction of the spin-liquid state for both small and 
large t±. In the latter region, the bonding and antibonding orbitals on a single rung separate 
to the extent that a state with two holons on a rung becomes higher in energy than two 
widely separated excitations each consisting of a hole in the bonding orbital of a rung. These 
quasiparticle excitations are described in RVB language as a bound holon-spinon pair, and 
were first observed by Tsunetsugu, et. al. [§J. The RVB picture suggests that the spin-liquid 
state is not really destroyed in this regime, only hidden by the spin degrees of freedom of the 
quasiparticles. This picture is quite compatible with the weak-coupling, band-separation 
picture mentioned above. 

In contrast, a band picture is not useful for small t±. In this regime, there is little 
difference in energy between the "staggered" and "resonant" valence bond configurations 
|16|| , and the average separation between two holons in a pair £/, becomes substantial. When 



£h becomes comparable to 1/5, where the doping 5=1 — (n), the pairs overlap, and the 
holons become essentially free. Since they are topological excitations residing on a single 
site, the free holons destroy the confinement mechanism responsible for the spin-liquid state. 
In this regime the system may resemble a doped system with an odd number of chains, where 
confinement never takes place. 

It is interesting to compare our results with several renormalization group calculations, 
which have also suggested the relevance of t± and the possibility of interchain pairing |8HTT 



Based upon a weak coupling renormalization group analysis, Fabrizio et al. [10] have sug- 
gested a possible U—t± phase diagram for the two chain Hubbard model. In the small t± 
region, they conjecture that there may be several phases in which pairing correlations are 



dominant. In one of these phases, called "SCI" in Ref. fllQfl , these pairing correlations have 
the same modified d x 2_ y 2 character we find. However, contrary to our results, they find 
a divergent spin-spin correlation function for q = (2kp, 0) rather than a spin gap state. 
Furthermore, they find this behavior for arbitrarily small values of t± and for generic, non- 
commensurate fillings. This differs from our findings that the pairing is associated with a 
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spin gap state and requires a finite value of t± > 0.5. Furthermore, we find that a doping of 
(n) near 1 is required for significant pairing. 



Recently, Khveshchenko and Rice [TTJ have carried out a renormalization group analysis 
near half-filling where umklapp processes become relevant. They argue that in the lightly 
doped system, the presence of the interchain hopping t± is sufficient to generate an anti- 
ferromagnetic exchange Jj_ leading to a spin gap, and that this state supports a modified 
d x 2^ y 2-\ike pairing. This scenario is certainly closer to what we find, except, as shown in 
Fig. |2|, we find that a finite value of t± is required to produce a spin gap. 
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Scalettar, and H. Tsunetsugu for useful discussions. R.M.N, and S.R.W. acknowledge sup- 
port from the Office of Naval Research under grant No. N00014-91-J-1143 and D.J.S. ac- 
knowledges support from the Department of Energy under grant DE-FG03-85ER45197 and 
the Program on Correlated Electrons at the Center for Materials Science at Los Alamos 
National Laboratory. The numerical calculations reported in this paper were performed at 
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FIGURES 

FIG. 1. (a) The charge gap Ac versus the inverse chain length L^ 1 for t± = 1, U = 8 and 
various fillings. The inset shows an enlarged view for (n) = 0.875. (b) The spin gap As versus L^ 1 
for t± = 1.0 at (n) = 1.0 and a number of t± values for (n) = 0.875. The lines are least squares 
fits to polynomials in L^ 1 . 

FIG. 2. The spin gap As versus t± for U = 8 and (n) = 0.875 calculated from the L — ► oo 
extrapolations of Fig. 1(b). The open circles show the magnitude of the pair field correlation 
function (AjAt) averaged over separations 8 < \i — j\ < 12 versus tj_ for the same parameters on 
a 2 x 32 lattice. 

FIG. 3. (a) The equal time rung-rung pair field correlation function (AjA^) versus \i — j\ 
plotted on a log-log scale for t± = 1.5 and for ladders of length L = 16, L = 24, and L = 32. The 
dashed line corresponds to the \i — j\~ 2 decay of the non-interacting pair field correlations, (b) 
The equal time pair field correlation function (AjAj) versus \i — j\ for various values of t±. 
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